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SUMMARY 


Flow field data for a double delta wing at low speed were used to determine 
the location of a vortex action point. The result was found to be consistent with 
what was determined for a delta wing. In supersonic flow, the action point location 
was determined empirically. For a wing with rounded leading edges, an assumption 
for initial vortex separation was shown to be equivalent to initial leading-edge 
bubble separation for airfoils. A theoretical formulation by the section analogy to 
determine the delayed vortex separation on a cambered wing with rounded leading edges 
was presented. The method of suction analogy was further shown to be applicable to 
predicting the body vortex lift. 


INTRODUCTION 


Since Polhamus introduced the method of suction analogy for plane delta wings 
with sharp edges in subsonic and supersonic flows (refs. 1 and 2), various 
extensions of the method have been proposed and used successfully. For wings with 
nonzero tip chord, Lamar introduced the effect of side-edge vortex lift (ref. 

3). When a vortex generated at the leading edge passes over some downstream 
planform area, additional lift is induced because of the vortex suction effect. 

This additional lift was incorporated into the method of suction analogy by Lamar 
as the augmented vortex lift (ref. 4). The latter may be positive as is the case 
for a cropped delta wing and a strake-wing configuration. Augmented vortex lift is 
negative for an arrow wing because of the lack of downstream area to receive the vortex 
action or to allow the reattachment of the vortex flow. For a noncambered wing, the 
vortex force was assumed to be acting at the leading edge and normal to the planform. 
However, for a positive-cambered wing, this assumption would produce too much 
thrust as compared with data. To allow for the determination of where the vortex 
force is acting, the concept of vortex action point was introduced (ref. 5). If 
the leading edge is rounded, the generation of vortex lift will be delayed. 

Kulfan presented a method to account for this effect for plane wings (refs. 6 and 
7). Another method for the effect of rounded leading edges was developed by 
Carlson by using available experimental data (ref. 8) . 

The method of suction analogy has also been employed to calculate the roll 
damping of slender wings (ref. 9). The extension of it to calculating all 
lateral-directional characteristics of slender wings was made in reference 10. In 
reference 10, the effect of vortex breakdown was also incorporated by using a 
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numerical scheme to extend the available data for delta wings. Finally, an 
unsteady lifting-surface theory was used in developing the method of unsteady 
suction analogy (ref. 11). The latter can be used for predicting dynamic 
stability derivatives of slender wings. 

In this paper, the concept of vortex action point is reconsidered for 

subsonic and supersonic flows. In addition, Kulfan's concept for plane wings with 
rounded leading edges is correlated with viscous flow calculations for airfoils. 
The extension of the method for cambered wings is described next. Finally, it is 
shown that the method of suction analogy can be applied to predict the vortex lift 
of slender bodies. 
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SYMBOLS 

local chord 

. 1 /2 
sectional leading-edge singularity parameter = limy(x)(x - x n ) 

Zb 

x+x, 

Zb 

drag coefficient 

sectional wave drag coefficient 

sectional lift coefficient 

lift coefficient 

pitching moment coefficient 

body normal force coefficient based on maximum cross-sectional area 
sectional normal force coefficient due to body vortex lift 
total augmented vortex normal force coefficient on a body 
pressure coefficient 

sectional leading-edge suction coefficient 
sectional side-force coefficient of a body 

total suction force coefficient produced by a body nose section 
sectional leading-edge thrust coefficient 

ratio of vortex-induced vertical velocity to the free stream (eq« (2)) 
vortex action point location on a body defined in eqs. (15) and (16) 

body nose length 
Mach number 
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r vortex action point location measured from the leading edge (eq, (4)), 

or body radius 

r Q leading-edge radius 

V velocity 

x s body axial station at which the vortex separation starts 

a angle of attack 

a 1 sectional angle of zero lift 

a g angle of attack of initial vortex separation 

0 sideslip angle or wave angle 

0_ oblique shock angle 

s 

Y ratio of specific heats or bounded vortex density 

A sweep angle 

p density 

0 angular 

0 A angular 

21 ) 

0 C angular 

Subscripts: 

£ lower surface 

Z e leading edge 

u upper surface 

°° free stream 


coordinate of a body (fig. 8) 

location of the augmented vortex action point on a body (eq. 
location of the body vortex action point (eq. (15) and fig. 10) 


CONCEPT OF VORTEX ACTION POINT 


In developing the concept of vortex action point, flow field data in cross- 
flow planes are needed, see figure 1(a). It is assumed that the velocity 
distribution around the vortex in the longitudinal plane parallel to the free 
stream is the same as that in the cross-flow plane (fig. 1(b)). If a control 

surface ( a) is taken through the vortex center as shown in figure 1(c), the vortex 
force should be equal to the vertical component of the force due to the momentum 
transfer through the control surface a. Since the vortex force acting on the wing 

is equal to c q c by the suction analogy, it follows that the force acting on the 
s 00 
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control volume is 


1 2 t + 

-c c — pV = J pV (V • dx) 
s 2 “ J z 


= /(PV. 


m 


pV ) V • dt 
Z out 


To find an average v, let 


( 1 ) 


I V I = V = fv (2) 

II z oo 

where f is a constant. It follows that for a unit span, 

/ ( pV - pV )V • dt = -fV [/pV dr + JpV dr] 

Z. Z 00 z . z 

in out m out 


= -fV 2p( fV )r (3) 


where the integrals have been replaced with the average values. From equations 
(1) and (3), it is found that 


r 



c c 
s 


(4) 


In reference 5, f was taken to be 0.5 so that r = c s c. To show that this choice 

of f is reasonable in subsonic flow, the flow field data in reference 12 for a 
double delta wing at station (1) are rearranged in figure 2. By numerically 
integrating the velocity distribution to evaluate the integrals in equation (3), 
an average velocity with f = 0.53 can be found. At station (2), the two leading- 
edge vortices have started to merge, so that the flow field data are not 
appropriate for the present purpose. 

In supersonic flow. Squire et al. demonstrated by a vapor screen method 
that the vortex region became much flatter than that in subsonic flow (ref. 13), 
as shown in figure 3. Unfortunately, flow field velocity data are not available 
in reference 13 and, in fact, have not been found so far for other wing 
planforms. Therefore, an empirical method t^ determine an appropriate f value was 
used. By assuming different values for 1 /4f and comparing the calculated C L and 
C D values with data for different planforms with leading-dege flaps, it was found 

that a value of 8.0 for 1/4f 2 produced the best results. In other words, f = 1/4/2. 
Therefore, unless additional data prove otherwise, the following relation for the 
vortex action point will be used in supersonic flow: 


r = 8(e g c) 


(5) 
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Some calculated results are compared with data in figures 4 and 5 (refs. 14 and 

15). It is seen that the theory predicts the results quite well. Note that the 
effect of rounded leading edges accounted for in figures 4 and 5 will be dis- 
cussed in the next section. The wave drag correction in figures 4 and 5 is to 
add the nonlinear effect to the predicted linear values based on the predicted 
sectional lift coefficients and the exact two-dimensional theory for a flat plate 
in supersonic flow. The wave drag correction is described in the Appendix. 


VORTEX LIFT ON CAMBERED WINGS WITH ROUNDED LEADING EDGES 


Kulfan assumes that on a slender wing the leading-edge vortex separation 
starts at an angle of attack at which the leading-edge drag equals the leading- 
edge thrust (refs. 6 and 7). To examine this assumption, experimental pressure 
data (ref. 16) and theoretical results from the Lockheed-Georgia airfoil code 
(ref. 17) were used. The airfoil negative pressure coefficient is integrated to 
give the leading-edge suction (c ), and the positive pressure coefficient is 

S £e 

integrated to produce the leading-edge drag (c ). The integration is over the 

d £e 

forward portion of the airfoil from the maximum thickness location if no 
separation bubble occurs. Otherwise, the integration is performed only over the 
region forward of the separation bubble. The results indicate that at the angle 
of attack (ct g ) at which the separation bubble first occurs, the leading-edge 

thrust is about equal to the leading-edge drag (ref. 18). For symmetrical 
airfoils, this is illustrated in figure 6. For cambered airfoils, some calculated 
results are presented in figure 7. As shown, the leading-edge suction at a g tends 

to be greater than the leading-edge drag. This is due to the pressure thrust 

generated on the forward camber. To remove this camber thrust in performing the 
pressure integration, the slope of the upper airfoil surface is reduced by the 
local camber slope and that of the lower surface is increased by it. The results 
indicated in figure 7 by rectangular symbols ( Q) show that Kulfan' s concept is 
still applicable for cambered airfoils if the camber thrust is removed from the 
calculated suction. In the thin airfoil or wing theory, this is always true 
because the calculated leading-edge thrust is concentrated at the leading edge and 
does not contain the camber thrust. 

Having established that Kulfan' s concept on the starting of a leading-edge 
vortex separation is related to the occurrence of leading-edge laminar separation, 
the next question is how this concept can be used in calculating the vortex lift 
on a cambered wing with rounded leading edges. The first task is to determine 
a g . For a cambered wing, the sectional leading-edge suction coefficient can be 
written as 


c = K( sina + a, ) 2 
s 1 


( 6 ) 


where K is a function of geometry and Mach number, and is the sectional angle 
of zero lift. In practice, can be found as follows. Let c g ' be the sectional 

leading-edge suction coefficient for the same wing without camber. Then 


69 



( 7 ) 


, • 2 
c ' = K sm a 
s 

From equations (6) and (7), a. can be found to be 

a, = [ (c /c ') 172 - i]sina (8) 

1 s s 

For a wing, the sectional leading-edge thrust coefficient (ref. 19) 
can be shown to be 


c = (tt/ 2)C 2 (1 - M 2 cos 2 A ) 172 /cosA (9) 

t 00 X© X© 

where C is the leading-edge singularity parameter and can be written as 


C = K' (sina + a ) 
1 


( 10 ) 


similar to equation (6). If C 1 is C at a g , then 


C = C(sina + a )/(sina + a ) 
1 si 1 


(11 ) 


The expression for the leading-edge drag can be found in reference 20. By 
equating the leading-edge drag to the leading-edge thrust at a g , it is obtained 
that 
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from which a g can be obtained: 


a 

s 


_.j sina + a 1 
sin [± 


(2 ^^"cosA^/d-M^cos^) 172 - Ul ] 


( 12 ) 


where r Q is the leading-edge radius. 

With a_, calculated, the remaining sectional thrust coefficient at a > a is 
s s 

then given by equation (9) with C replaced by C2 , where 

C = C[sin( a - a ) + a ] /(sina + a ) (13) 

2 si' 1 
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Note that c. must be converted to the sectional suction coefficient (c„ ) before 

the latter is assumed to become the vortex lift through the method of suction 
analogy. The relation between c. and c g for a cambered wing was derived in 

reference 5. 

The above consideration has been applied to wings in both subsonic and 
supersonic flows with good success. In the supersonic flow, the wing must have a 
subsonic leading edge to produce the vortex lift. Some correlation with data was 
shown earlier in figures 4 and 5. 


BODY VORTEX LIFT 


If the aspect ratio of a slender thick wing is reduced, eventually it becomes 
a slender body. If the method of suction analogy is applicable to the former, it 
should also be applicable to the latter. Based on this understanding, the 
following method for calculating body vortex lift was developed. The method is 
based on the following assumptions and procedures. 


(a) The attached-flow solution is obtained with the axial distribution of G. 
N. Ward's vortex multiplets (ref. 21). The boundary condition is 
satisfied on the body surface. 

(b) At any axial station, vortex separation starts at a circumferential 
location where C is minimum and negative. This assumption has been 

p 

shown to be reasonable (ref. 22). At low angles of attack, C may be 

Jr 

positive everywhere near the nose. In this case, no vortex separation 
is assumed to occur. In reference 22, the axial station (x s ) at which 

the vortex separation starts must be assumed or given by experimental 
data. Examples of calculated pressure distribution with and without 
vortex separation are illustrated in figure 8. 

(c) At any axial station, the side-force component of the negative in the 

region assumed to have vortex separation is integrated to produce a 
sectional side-force coefficient (see fig. 9). 

6 . 

mm p 

c = — ; — - / r(x)C sin0d0 (14) 

sf r(x) 1 p 

The side force obtained is assumed to be the suction force produced by 
the separated vortex. This suction force is assumed to be acting at 0 C 

(fig. 10) where 


0=0. - A0 (15) 

c mm p 


A0 = h c 

sf 


( 16 ) 
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Equations (15)- (16) imply that the vortex action point is located at a 
distance from 0 m £ n being proportional to the suction force: 


As = rA0 = h r c 

sf 

Based on equation (4), h should be 


(17) 


h = 1 /4f 2 


(18) 


To determine an appropriate value for f, data in reference 23 as 
presented in figure 11 are used. Applying the same method as used in 
figure 2, f is determined to be 0.70. For simplicity, f will be taken 
to be l/(2) 1 / 2 in subsonic flow. In supersonic flow, no data were 
available for correlation so that an empirical value will be used. In 
summary, the following values for h will be used in the present method: 


0 . 

,5 

, M 

< 

1.0 



00 



3. 

,0 


> 

1.0 


oo 


With 0 C calculated from equation (15), the sectional normal force 
coefficient due to the vortex suction is given by 


c = c COS0 
n,v sf c 


( 20 ) 


(d) Similar to a wing, the augmented body vortex lift exists whenever the 
planview of a body is not of the delta type. The concept used in 
calculating the augmented vortex lift for a wing (ref. 4) is also 
applicable for a body. Thus, if C g £ N is the total suction force 

coefficient from the nose portion, then the augmented normal force 
coefficient (C N vA ) is given by 


C = c C _ cosG /£ 
N,vA sf,N A N 


(21 ) 


where c is the length over which the nose vortex passes and is the 

nose length. The term 0 ft is the location of augmented vortex action 
point and is assumed to be equal to 0 c at the body shoulder. 


Calculated results for an ogive cylinder at = 0.3 are compared 

with data in figure 12, The effect of Reynolds number is seen to increase 
the loading. In references 22 (ref. 24), different axial locations of 



initial separation were assumed for different Reynolds numbers. The 
present results are seen to agree well with the data of higher Reynolds 
number. This is expected because an inviscid theory, such as the 
present one, is to simulate the flow field of infinite Reynolds number. 

In figure 13 (ref. 25), results for a circular-arc-cylinder body at 

Moo = 1 .6 are presented. It is seen that the calculated results agree well 
with data up to a = 32 deg. Above a = 32 deg at M OT = or 

a > 15 deg at M =2.3 (not shown), the normal force is always 
00 

underpredicted. This is probably because at high M and/or high angles 

■ ' ■ - . 00 

of attack, the present linear theory cannot predict accurately the upper 
surface expansion and strong shock effect on the lower surface. 


CONCLUDING REMARKS 


The method of suction analogy, originally developed for a plane wing, was 
shown to be applicable to cambered wings by using the concept of vortex action 
point. For a wing with rounded leading edges, the method is still applicable if 
the delay in initial vortex separation is accounted for. The latter can be 
calculated by Kulfan's method. Kulfan's method was shown to be related to the 
leading-edge laminar separation. Extension of the suction analogy to predicting 
the body vortex lift has also been presented. 
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APPENDIX 


CORRECTION FOR NONLINEAR WAVE DRAG 


The correction for nonlinear wave drag is based on the difference in 
predicted results for a two-dimensional flat plate by the linear and exact 
theories. The linear theory shows that the drag coefficient is given by 

c, = 4cosa sin 2 a/(M 2 - l) 1 ^ 2 (A.1) 

d 00 

w 

where the a- terms are not linearized, so that the linearization is with respect to 
the compressibility effect only. For a three-dimensional wing, a in equation 
(A.1) is the one associated with the predicted sectional lift coefficient (c^) as 
follows: 


sin 1 [c „ (M 

X. 0 


- 1) 1/2 ]/4 


(A. 2) 


The exact theory for a flat plate in supersonic flow shows that the upper 
surface pressure is given by the Prandtl-Meyer solution (ref. 26, p. 383): 


a .. 3l l?„ ) 1f/(Y ' 1l /( 


2 

sin 3 


P „ 2 y - cos 2 3 

u ym 


Y - cos 2 3 


fL_jY/(Y-1) _ u 


(A. 3) 


where 3 = sin 1 ( 1 /M ) and 3 is found from 


„ ,Y + 1,1/2 -1 r ,Y + 1,1/2 

-a + 3 - ( ^ _ — ) tan [ ( - — ■ ^ ) tan33 


■ e. - <74T ),/2t “' 1[| 7TT |1/2tanB - 1 


( A. 4 ) 


The lower surface pressure is given by the oblique shock solution (ref. 27, p. 
86-88): 


C = 2sin3 sina/cos(3 - a) (A. 5) 

P £ 

where 3 S is obtained from 

2M 2 sin 2 3 - 2 

oo s 

tana = cot3 r (A. 6) 

S (Y+1)M‘ c '- 2M sin g + 2 

oo oo g 

From equations (A. 3) and (A. 5), the pressure differential and the wave drag become 
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(A, 7) 


AC 


P 



C 


c = AC sina (A. 8 ) 

d p 

w 

The difference between the values given by equations (A. 8) and (A.1) 
represents the correction to be added to the predicted sectional drag coefficient 
by the linear lifting surface theory. 
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VORTEX 

REGION 


Figure 3 


igure 4. 



. Vortex flow regions on a 65-deg 
and supersonic speeds. a = 8.1 


plane wing at subsonic 
deg and R g = 1.5 X 10®. 


THEORY - ROUNDED EDGES 




Longitudinal characteristics of a delta wing of aspect ratio 
4/3 at M =1.8. 
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Figure 7. Camber effect on leading-edge suction at angles of attack at 


which separation bubble first occurs. M =0.17, 


R = 5.8 X 10 
e 



e, deg 


Figure 8. Calculated pressure distribution on an ogive cylinder at 
a = 10 degrees and M =0.3. 
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Figure 12 . 
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11. Flow field data from ref. 23 on a tangent ogive body at 
x/D = 2.6, a = 45 deg and = 1.5 X 10^. 


PRESENT THEORY 

VORTEX FLOW 

ATTACHED FLOW 



Normal force and pitching moment coefficients for an ogive-cylinder 
body at M =0.3. Moment center at body nose. Reference area = 

“ oo 

base area and reference length = d. 








PRESENT THEORY 

VORTEX FLOW 

ATTACHED FLOW 



Figure 13. Normal force and pitching moment coefficients for a circular-arc- 
cylincer body at = 1.6. Moment center at body nose. Reference 
area = base area and reference length = body length. 


84 



